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Abstract. - We report an inelastic neutron scattering (INS) study under a magnetic field on 
the frustrated molecular spin cluster V15. Several field-dependent transitions are observed and 
provide a comprehensive understanding of the low-energy quantum spin states. The energy 
gap 2Ao « 27(3)/xeV between the two lowest S — 1/2 Kramers doublets is unambiguously 
attributed to a symmetry lowering of the cluster. The INS data are mapped onto an S=l/2 
Antiferromagnetic Heisenberg triangle with scalene distortion. A quantitative description of 
the wavefunction mixing within the ground state is derived. 



Magnetic frustration operates when all the bonds in a magnetically coupled system can- 
not be satisfied simultaneously. The ground state is then best described by a superposition 
of quantum states with well defined probability of occurrence. However, small irregularities 
like lattice strains, structural disorder or quantum fluctuations can relieve part or all of the 
frustration, leading to a stabilised ground state with a gap to less favourable ground state 
configurations (order by disorder principle [1]). The field of frustrated magnetism is extremely 
active, encompassing low-dimensional materials [2,3], spin glasses or pyrochlores and Kagome 
lattices [4] . Most of these materials are extended systems exhibiting short-range or quasi long- 
range order, a situation rarely encountered in molecular magnets (or spin clusters) where each 
such cluster in the lattice is magnetically well isolated from its neighbours due to the presence 
of surrounding ligands. This magnetic shielding allows the study of the individual behaviour 
of a finite number of interacting magnetic ions. Beside being ideal candidates to study funda- 
mental processes like quantum tunnelling or quantum coherence at the nanoscale [5] , molecular 
magnetic clusters, with their well defined nuclearities and topologies, also constitute quantum 
systems in which geometrical frustration effects at the molecular level can be addressed. 
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A prominent example of a magnetically frustrated cluster is given by the the polyoxo- 
vanadate complex K 6 [V 1 7 5 y As 6 042(-D20)] • 8D 2 (hereafter V15) which contains 15 V 4+ spins 
(S — 1/2) distributed over two hexagons capping one triangle with global spherical shape [6-8] . 
The clusters have D3 symmetry without consideration of the hydrogen positions (inside and 
outside the cluster) and make up a molecular crystal with trigonal symmetry. The V A+ ions 
are antiferromagnetically (AFM) coupled to their neighbours via oxo-bridges. Within the 
hexagons, the AFM couplings are very strong (10-20 meV [9]) and the spins on the triangle 
are coupled to the spins of the hexagons via frustrated exchange couplings, but there are no 
significant direct exchange pathways between the triangle spins. The coupling between the 
triangle spins occurs indirectly through the hexagon spins. At low temperatures, the spins 
on the hexagons are quenched in a singlet (S=0) state, and we are left with the three tri- 
angle spins. The ground state is then made of two S=l/2 Kramers doublets separated from 
the S=3/2 quartet state by about 0.315 meV (3.7K) [10, 11]. Inelastic neutron scattering 
(INS) has shown that the two S=l/2 Kramers doublets are split by a gap 2Ao ~ 35^eV [10], 
and low-temperature magnetisation data were analysed in terms of a two-level Landau-Zener 
model with a gap of about 7 — 8/ieV (80 — 100 mK) [11]. However, the microscopic origin of 
these gaps remained unclear. Recent theoretical work suggested that Dzyaloshinskii-Moriya 
(DM) interactions might be responsible for the gap opening [12]. 

In the present Letter, we report an INS study under a magnetic field of V15 and show that 
deviations from trigonal symmetry are responsible for the observed phenomena, not DM inter- 
actions. An analysis of both the energy levels and the Q-dependence enables us to characterise 
the distortion and to address the nature of the ground state. 

The INS experiment was performed on the recently upgraded time-of-flight spectrometer 
IN5 at the Institut Laue-Langevin (ILL, France) using cold neutrons of wavelengths A = 7.5, 9 
and 11 A. Data were collected at temperatures between 40mK and 50mK and corrected for 
the background and detector efficiency. The magnetic field was provided by a 2.5 Tesla 
superconducting coil. The instrumental resolution at the elastic line (r, Full- Width at Half- 
Maximum) is T = 28/ieV at 7.5A, Y = 18 - 21^eV at 9.0A and V = 12/ieV at llA. We 
used a 4.6 g polycrystalline powder sample of fully deuterated V15 placed under Helium in a 
rectangular flat Aluminum slab. 

Figure^shows INS spectra obtained at 7.5A and 9.0A for different values of the magnetic 
field. At H=0T, only one transition can be observed at an energy of w 0.335 meV with 
r w 41/xeV. This width is 1.5 broader than the instrumental resolution and is intrinsic [14]. 
As the field is switched on, satellite peaks appear symmetrically on each side of the main peak 
but with different intensities. At IT, there is new intensity at about 0.12 meV. To better 
characterise it, the difference between the H=0.5/1T and 0T data is shown in fig. ^ along 
with a higher resolution (A = 11 A) spectrum obtained at 0.65T. There are two peaks separated 
by about ss 25 — 30/ieV, and their energies increase linearly with magnetic field. Their widths 
are found to be between 17 and 19/ieV, in the range of the instrumental resolution, in contrast 
to the higher-energy peaks shown in fig. ^ Compiling the information derived from fig. ^ 
it is possible to draw the field-dependence of the INS transitions (see fig. Et) where the five 
transitions are labeled (I) to (V). 

To model our data, we consider the S = 1/2 Heisenberg Antiferromagnetic (HAFM) model 
on a triangle [7, 8]: 

7~Lo = J12S1S2 + J23S2S3 + J13S1S3 + Hh , (1) 

where Si, S2, S3 denote the spin operators 1,2 and 3, respectively and is the Heisenberg 
exchange parameter between spins i and j. Jin = t L BH{giS\ + 92 S2 + .93 S3) is the Zeeman 
interaction where the (^-tensors are assumed to be diagonal. The Zeeman splitting of the 




Fig. 1 - (a) INS spectra at 0,0.5 and 0.75T (A = 7.5A). (b) INS spectra at and LOT (A = 9.0A) 
(c) INS spectra at 0.65T (A = 11. OA) and difference plots between non-zero magnetic field (H=0.5T, 
LOT) and zero-field data obtained at A = 9.0A. Data in the Q-range between 0.2A -1 and 1.1-1.4A -1 
were grouped together to improve statistics. Solid lines are best fits to the data using Gaussian line 
shapes and a background. Peaks are labeled as discussed in the text. 



energy levels does not depend on the relative orientation to the trigonal axis for isotropic 
couplings. For the equilateral HAFM triangle case (Jy — Jo > 0), the ground state consists 
of two degenerate S=l/2 Kramers doublets separated from the S = 3/2 excited state by 
3Jo/2. The wave functions of the two degenerate S=l/2 Kramers doublets are given by 

*o 2 = |0) i> -^-h ^ ano - ^i 2 = I-*-' 3>^-3 ^ m ^ ne basis \S\2, S, M > or any linear combination 
of them. If we assume inequivalent couplings, there is a gap 2Ao between the two Kramers 

doublets but no splitting in the S — 3/2 state. For instance, in the isosceles case (J12 = J, 

-l-i 

Jl3 = J23 = J') with J > J', the energy gap becomes 2A = J — J' and ^ 2 is the lowest 
doublet. 

From the field-dependence of the INS transitions shown in fig. we can immediately 
construct the energy diagram shown in fig. |2h corresponding to a distorted S=l/2 HAFM 
triangle with two S=l/2 doublet states and one S=3/2 quartet state. Note that at T = 45 
mK, only the lowest doublet is populated in zero- field. Transitions (I) and (II) are intra- 
doublet transitions with energies hu>i = h = g^sH, Tuon = h + 2Aq, and transitions (III) to 
(V) are transitions from the lowest doublet to the S=3/2 sublevels with faom = 3Jo/2+Ao— h, 
Ticujv — 3Jo/2 + Ao and fitoy = 3 Jo/2 + Ao + h. The transition two — 2Aq between the two 
doublets at zero-field could not be observed due to the width of the elastic line. The solid 
lines in fig. UK corresponds to best fits with g = 1.98(3), 2A = 27(3)^eV » 0.31(4)K and 
Jo = 212(2)/ieV « 2.46K. The g-factor value is in good agreement with EPR [7, 13] and 
millimeter range spectroscopy [15], the gap 2Ao is in relatively good agreement with previous 
but less accurate INS measurements [10], and the value of Jo is very close to the one inferred 
from magnetisation [11] and INS [10]. So far, the isosceles triangle scenario successfully 
explains the existing data, but there is no reason, a priori, to choose an isosceles triangle: 
Lattice distortions, being static or dynamic, will most probably generate a scalene triangle 
situation. From inspection of the energy levels it is not possible to discriminate between 
isosceles and scalene triangles. We demonstrate now that the INS intensities can discriminate 
between the two models. 
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Fig. 2 - (a) Field-dependence energies of the observed INS peaks (I) to (V). Solid lines are linear 
fits to the data as discussed in the text, (b) Zeeman splittings and assignment of the observed INS 
transitions. 



The differential magnetic cross-section for a transition between initial state \S\ 2 ,S,M > 
with energy E s and final state \S[ 2 , S', M' > with energy Ef is given by [16]: 

x (S 12 ,S,M\S?\S[ 2 ,S',M')(S' 12 ,S',M'\S?\S 12 ,S,M) x6(hu + E s - E f ) , 



with 

B Q = ^P-^.F\Q).e-^. (3) 

N is the number of magnetic centres in the sample, k and k' are the initial and final neutron 
wave-vectors and Q = k — k' is the scattering vector, exp[— 2W{Q)\ is the Debye- Waller 
factor, F(Q) is the magnetic form factor of the V 4+ ions, Hi is the position of the ith V 4+ 
ion in the triangle, a = x, y or z, hu> is the energy transfer and Z the partition function. The 
matrix elements (S'12, S, M\Sf \S[ 2 , S', M') are evaluated using irreducible tensor operator 
(ITO) methods [17]. Only the triangular model is considered. In the equilateral case, the 
zero- field and powder averaged cross-section for the transition | S'12, S > to | S' 12 ,S' > is [18]: 

,Pa sJ^^m^M^S^S'), (4) 



dQ.dE v V Q R 

where R is the V 4+ -V 4+ separation in the triangle and M is a number that depends on the 
initial and final state quantum numbers. One can show that it makes no difference to the 
INS intensities in zero magnetic field whether we have S'12 = or S12 = 1 in the lowest 
doublet. In order to compare the theoretical predictions with our experimental data, we 
need to consider the effect of the magnetic field. This is done by using the Wigner-Eckart 
theorem [18]. This leads to INS intensities depending explicitly on M, M' with selection rules 
M — M' = 0, ±1. We obtain the following intensity ratios for the doublet-quartet transitions 
(III:IV:V) = (3:2:1). The intra-doublet transition with AS 12 = S 12 - S' 12 = ±1 has the same 
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Fig. 3 - (a) Q-dependence of the transitions III + IV + V measured at zero-field, shown in fig. 
(b) Q-dependence of transition (I) at IT and A = 9A (open squares) along with the III+IV+V 
transitions measured at 0T and at the same wavelength. Lines (dashed and solid) are best fits using 
the calculated Q-dependence (eq. and the V — V distance R — 6.92A. Dotted lines correspond to 
theoretical curves with a 2 — 1 and b 2 — 1 (see text). 



relative intensity as transition (IV), whereas the AS12 = transition is calculated to be three 
times more intense. Comparison with the data in fig.^shows good, but not perfect agreement 
for the doublet-quartet transitions, as (III) and (IV) have almost equal intensity, instead of 
the calculated ratio 3/2. From the relative intensity of the intra-doublet transitions we can 
tentatively assign transition (I) to have AS12 = and transition (II) to have AS12 = ±1. We 
already note at this point, however, that transition (I) is weaker than the sum of (III), (IV) 
and (V). 

To get further insight into the nature of the transitions, in particular the intra-doublet 
ones, we now consider their Q-dependence. Fig. shows the Q-dependence of the sum of 
(III), (IV) and (V), measured at zero- field. The agreement with the Q-dependence calculated 
with Eq.Q]and the V-V distance R = 6.92A is very good, a confirmation that the S=3/2 state 
is essentially unperturbed. Fig. [SJ^ shows the same transitions, i.e. (III+IV+V) measured at 
zero- field and A = 9A (full circles), in comparison with the Q-dependence of transition (I) at 
IT at the same wavelength (open squares) [19]. While the doublet to quartet transition shows 
the same Q-dependence as in fig. the intensity of (I) is much less Q-dependent, almost 
flat. Theoretically, an intra-doublet transition with AS12 = ±1 would be modulated by the 
[1 — sin(Qi?)/Qi?)] factor in eq.0J On the other hand, a AS12 = transition does not not have 
this factor and its Q-dependence is essentially flat. We have already tentatively assigned band 
(I) to a Zeeman transition within the same doublet (AS 12 — 0), and the Q-dependence now 
confirms this. However, we also note, and fig. Ob makes it very clear, that the overall intensity 
of peak (I) is significantly smaller than the sum of (III+IV+V). The pure AS12 = transition 
would have the same intensity, and the pure AS12 = ±1 would be three times weaker. This 
clearly suggests that the two doublets Vf/j^ and are mixed in V15 as expected for a scalene 
triangle. Defining the lower-lying doublet as fl^ = aVl/^ + b^f (with a 2 + b 2 = 1), we get the 
following Q-dependent intensity for transition (I): 

Ii(Q)=hF 2 (Q) 

where I is an intensity factor, kept fixed to the value determined from the intensity of 
(III+IV+V) at the same wavelength. The dotted lines in fig. |3b represent the disentangled 
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situations ([a, 6]=[1,0] or [0,1]) and the full line is a best fit to the data corresponding to 
a 2 = 0.4 and b 2 = 0.6. The mixing of the S12 = and S12 = 1 states is quite substantial. One 
parameter set that produces this situation is J12 = 0.21 meV, J23 = 0.23 meV and J13 = 0.20 
meV. This set is not unique as we only have access to a 2 and b 2 , but it gives a clear idea of 
the exchange coupling variation caused by the distortion. 

We now briefly discuss the relevance of Dzyaloshinskii-Moriya (DM) interactions [12]. The 
effect of DM interactions in the S=l/2 AFM triangle has been treated in some detail [20] and 
we give here the main results. The DM interaction introduces off-diagonal matrix elements 
which lead to a first order splitting 2A of the doublets. If we assume perfect triangular 
symmetry, we have 2A = d z /^/3 [21] where d z is the DM parameter. If the splitting in V15 was 
due to DM interactions this would lead to d z = 47/icV. In contrast to the pure Heisenberg case, 
the energy levels depend explicitly on the relative orientation of the clusters to the magnetic 
field in this case [20]. In a powder measurement, the overall spectrum will reflect this field 
dependence by a broadening of all the peaks, in particular the intra-doublet transitions. This 
is clearly not what we observe experimentally. In addition, the theoretical energy difference 
2A in the DM model between transitions (I) and (II) is field dependent, decreasing from 
27/ieV at 0T to 17/ieV at IT for a d z value of 47/ieV. Again, this is completely incompatible 
with our experimental results, which show a field independent energy difference between (I) 
and (II) as well as a linear field dependence of transitions (III), (IV) and (V), see fig. [21 DM 
interactions can thus definitely be ruled out as the origin of the ground state splitting in V15. 

The present study demonstrates that the low-energy properties of V15 are accurately de- 
scribed by a triangle model with scalene distortion. In particular, the mechanism generating 
a gap 2Ao = 27(3)/xeV between the two Kramers doublets is unambiguously established us- 
ing both energy and wavefunction information provided by INS. The symmetry lowering is 
already established from the field dependence of the energy levels. However, only inspection 
of the intensities and their Q-dependence leads to a detailed knowledge of the wavefunction 
mixing within the ground state resulting from the scalene distortion of the triangle. An ob- 
vious source of symmetry lowering is provided by the water molecule located in the center of 
the spherical cavity of V15. Another one is given by some disorder (partial occupancy) on 
the water structure of the lattice [7,8,22]. This study shows that the impact of very small 
structural perturbations to relieve magnetic frustration can be sizeable. V15 constitutes a 
example of a nanometer-scale system with maximised magnetic frustration (triangle) where 
order, i. e. a stabilised ground state with a gap to less favourable " ground states" , is induced 
by a small structural distortion (disorder). 
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